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Abstract 

This document outlines a set of simplified models for dark matter and its interactions with Stan¬ 
dard Model particles. It is intended to summarize the main characteristics that these simplified 
models have when applied to dark matter searches at the LHC, and to provide a number of useful 
expressions for reference. The list of models includes both s-channel and t-channel scenarios. For 
s-channel, spin-0 and spin-1 mediation is discussed, and also realizations where the Higgs particle 
provides a portal between the dark and visible sectors. The guiding principles underpinning the 
proposed simplified models are spelled out, and some suggestions for implementation are presented. 
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1. Introduction 


Gravitational effects on astrophysical scales give convincing evidence for the presence of dark 
matter (DM) in Nature, an observation that is strongly supported by the large-scale structure of the 
Universe and measurements of the cosmic microwave background [T]. While the existence of DM 
thus seems well established, very little is known about the properties of the DM particle(s). To shed 
light on this question, three classes of search strategies are being employed: (i) direct detection 
in shielded underground detectors; {ii) indirect detection with satellites, balloons, and ground- 
based telescopes looking for signals of DM annihilation; (Hi) particle colliders aiming at direct DM 
production. Despite this intense effort, DM has so far proven elusive. In the coming years, direct 
and indirect detection will reach new levels of sensitivity, and the LHC will be operating at 13 TeV 
centre-of-mass energy after a very successful 8 TeV run. These upcoming experiments will provide 
crucial tests of our ideas about DM, and have great potential to revolutionize our understanding 
of its nature. 

Dedicated searches for DM candidates represent an integral part of the physics programme at 
the LHC. The minimal experimental signature of DM production at a hadron collider consists of 
an excess of events with a single final-state object X recoiling against large amounts of missing 
transverse momentum or energy (^t)- In Run I of the LHC, the ATLAS and CMS collaborations 
have examined a variety of such “mono-A” signatures involving jets of hadrons, gauge bosons, top 
and bottom quarks as well as the Higgs boson in the final state. A second class of signatures 
that has been studied in depth arises from the production of “partner” particles that decay to DM 
and Standard Model (SM) particles, which usually leads to rather complex final states (for a review 
of the experimental status after LHC Run I, see for instance [2]). 

In order to interpret the cross section limits obtained from the LHC $x searches, and to relate 
these bounds to the constraints that derive from direct and indirect detection, one needs a theory 
of DM. In fact, as illustrated in Figure one can construct not just one, but a large number of 
qualitatively different DM models. Collectively these models populate the “theory space” of all 
possible realizations of physics beyond the SM with a particle that is a viable DM candidate. The 
members of this theory space fall into three distinct classes: 

(I) On the simple end of the spectrum, we have theories where the DM may be the only accessible 
state to our experiments. In such a case, effective field theory (EFT) allows us to describe the 
DM-SM interactions mediated by all kinematically inaccessible particles in a universal way. 
The DM-EFT approach [3101 El [6l [71 [HI [9l [10] has proven to be very useful in the analysis of 
LHC Run I data, because it allows to derive stringent bounds on the “new-physics” scale A 
that suppresses the higher-dimensional operators. Since for each operator a single parameter 
encodes the information on all the heavy states of the dark sector, comparing LHC bounds 
to the limits following from direct and indirect DM searches is straightforward in the context 
of DM-EFTs. 

(H) The large energies accessible at the LHC call into question the momentum expansion under¬ 
lying the EFT approximation [giini Emni [13101151 diiiiTi , and we can expand our level of 
detail toward simplified DM models (for early proposals see for example [THUTOl 12011^1221123] L 
Such models are characterized by the most important state mediating the DM particle in¬ 
teractions with the SM, as well as the DM particle itself. Unlike the DM-EFTs, simplified 
models are able to describe correctly the full kinematics of DM production at the LHC, be¬ 
cause they resolve the EFT contact interactions into single-particle s-channel or t-channel 
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Figure 1: Artistic view of the DM theory space. See text for detailed explanations. 


exchanges. This comes with the price that they typically involve not just one, but a handful 
of parameters that characterize the dark sector and its coupling to the visible sector. 

(Ill) While simplified models capture some set of signals accurately at LHC energies (and beyond), 
they are likely to miss important correlations between observables. Complete DM models 
close this gap by adding more particles to the SM, most of which are not suitable DM 
candidates. The classical example is the Minimal Supersymmetric SM (MSSM), in which 
each SM particle gets its own superpartner and the DM candidate, the neutralino, is a weakly 
interacting massive particle. Reasonable phenomenological models in this class have of order 
20 parameters, leading to varied visions of DM. At the same time, they build-in correlations 
from symmetry-enforcing relations among couplings, that would look like random accidents 
in a simplified model description. Complete DM models can in principle answer any question 
satisfactorily, but one might worry that their structure is so rich that it is impossible to 
determine unambiguously the underlying new dynamics from a finite amount of data (“inverse 
problem”) [2lj. 

Given our ignorance of the portal(s) between the dark sector and the SM, it is important that 
we explore all possibilities that the DM theory space has to offer. While the three frameworks 
discussed above have their own pros and cons, they are all well-motivated, interesting, and each 
could, on its own, very well lead to breakthroughs in our understanding of DM. Ignoring whole 
“continents” of the DM theory landscape at Run II, say EFTs, would be shortsighted, and might 
well make it impossible to exploit the full LHC potential as a DM discovery machine. 

In recent years, a lot of progress has been made in exploring and understanding both DM-EFTs 
and a variety of complete models. The same cannot (yet) be said about simplihed models that 
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bridge between the two ends of the spectrum in theory space. Following the spirit of |25| I26|. we 
attempt in this document to lay the theoretical groundwork that should be useful for the DM@LHC 
practitioner. We begin in Section by discussing the general criteria that a simplified DM model 
should fulfill to make it useful at the LHC. This section contains in addition an explanation of the 
concept of Minimal Flavor Violation (MFV) |27l [28l [29l [30] and its importance to model building 
as well as a brief note on the relevance of the spin of the DM particle for LHC searches. Simplified 
spin-0 s-channel models are then described in Section]^ Since these scenarios can be understood 
as limiting cases of Higgs portal models, we provide in Section]^ a summary of the most important 
representatives of these theories. Section]^ is devoted to simplified spin-1 s-channel models, while 
Section deals with t-channel scenarios. To make the work self-contained, we not only discuss the 
LHC phenomenology of each simplified model, but also provide the relevant formulae to analyze 
the constraints from direct detection and annihilation of DM. We conclude and provide an ontlook 
in Section [3 

2. Criteria for Simplified Models 

For a simplified DM model to be useful at the LHC, it should fulfill the following three criteria: 
(i) it should be simple enough to form a credible unit within a more complicated model; {ii) it 
should be complete enough to be able to describe accurately the relevant physics phenomena at 
the energies that can be probed at the LHC; {iii) by construction it should satisfy all non high-p-r 
constraints in most of its parameter space. 

One way to guarantee that these three criteria are met consists in putting the following require¬ 
ments/restrictions on the particle content and the interactions of the simplified model: 

(I) Besides the SM, the model should contain a DM candidate that is either absolutely stable 
or lives long enough to escape the LHC detectors, as well as a mediator that couples the 
two sectors. The dark sector can be richer, but the additional states should be somewhat 
decoupled. A typical mass spectrum is sketched on the left in Figure]^ 

(H) The Lagrangian should contain (in principle) all terms that are renormalizable and consistent 
with Lorentz invariance, the SM gauge symmetries, and DM stability. However, it may 
be permissible to neglect interactions or to study cases where couplings are set equal to 
one another. If such simplifications are made, one should however try to verify that these 
approximations do not result in a very different DM phenomenology and they should be 
spelled out clearly in the text and on all relevant plots. 

(HI) The additional interactions should not violate the exact and approximate accidental global 
symmetries of the SM. This means that the interactions between the visible and the dark 
sector should be such that baryon and lepton number is conserved and that the custodial and 
flavor symmetries of the SM are not strongly broken. 

Simplified models are thus specifically designed to involve only a few new particles and interac¬ 
tions, and many of them can be understood as a limit of a more general new-physics scenario, where 
all but the lightest dark-sector states are integrated out. By construction, the physics of simplified 
models can therefore be characterized in terms of a small number of parameters such as particle 
masses and couplings. While simplified models are clearly not model-independent, they do avoid 
some pitfalls of DM-EFTs. In particular, they allow one to correctly describe the kinematics of DM 
production at the LHC, by virtue of the dynamical mediator(s) that they contain. Conversely, by 

5 


lOTeV -I- 

ITeV 
100 GeV 




t,h,Z,W 


X 


lOTeV + 

ITeV 
100 GeV 




^3 


t,h,Z,W 


X 


Figure 2: Left: Schematic mass spectrum of a simplified DM model. In the case considered, the DM particle x is 
lighter than the heaviest SM particles t, h, Z, W. The lightest mediator state is called Zi and can be produced on- 
shell at the LHC. The remaining dark-sector states Z 2 and Z 3 are separated by a mass gap from Zi and inaccessible. 
Right: The EFT limit of the simplified model with a decoupled mediator Zi. See text for further details. 


making the mediator(s) sufficiently heavy the EFT framework can be recovered. The latter feature 
is illustrated on the left-hand side of Figure 

2.1. Note about Flavor and CP Violation 

The requirement (III) deserves further explanations. The SM posseses both exact and approx¬ 
imate global accidental symmetries. The former (baryon and lepton number) are conserved at the 
renormalizable level, while the latter (custodial and flavor symmetries) are broken by quantum 
effects, but parametrically small in the sense that they become exact global symmetries when a 
parameter or a number of parameters are set to zero. New physics will generically not respect these 
accidental symmetries and, as a result, its parameter space will be severely constrained: the new 
interactions are required to be weak or the new states have to be heavy (or both). 

A systematic way to curb the size of dangerous flavor-violating and CP-violating effects consists 
in imposing MFV. Loosely speaking the idea behind MFV is that the general structure of flavor¬ 
changing neutral current (FCNC) processes present in the SM is preserved by new physics. In 
particular, all flavor-violating and CP-violating transitions are governed by the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix. This basic idea can be formalized and formulated in an EFT [30]. Employ¬ 
ing the EFT language, a new-physics model satisfies the MFV criterion if the additional interactions 
in the quark sector are either invariant under the global SM flavor group Gq = U{3)g x U (3)^ x U (3)^, 
or any breaking is associated with the quark Yukawa matrices Y“ and The notion of MFV can 
be also be extended to the case of CP violation and to the lepton sector — although for leptons its 
definition is not unique, if one wants to accommodate neutrino masses. 


2.1.1. MFV Spin-0 s-Channel Models 

To understand which restrictions MFV imposes on the flavor structure of simplified models, 
we work out some examples relevant for the discussions in later sections. We begin with a very 
simple model in which DM is a real scalar (gauge and flavor) singlet x and the SM Higgs doublet 
H provides a portal to the dark sector of the form x^l-Hp (the most important phenomenological 
implications of this scenario will be discussed in Section 4.2). Following the notion of MFV, the 
interaction terms between the mediator and the quark fields should be either invariant under Gq 
or break it only via Y“ or Y'^. Given the transformation properties q (3,1,1), u (1, 3,1) and 
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d ~ (1,1,3), it follows that the combination qu of left-handed and right-handed quark fields breaks 
U{3)q X U{3)u, while the bilinear qd breaks U{3)q x U{3)d- This means that we have to go with 
the second option. In terms of gauge eigenstates, we write 

T D - ^ (^{Y^)^jqiHuq + {Y%qiHdq + h.c.) , ( 1 ) 

where i,j runs over the three quark families, Ha = ^abH^ with a,b = 1,2 and the two terms involve 
the Higgs fields to make them SU{2)l x U{1)y gauge invariant. Notice that the above interactions 
are invariant under Gq, if the Yukawa matrices are promoted to non-dynamical helds (spurious) 
with the following transformation properties Y“ ~ (3,3,1) and ~ (3,1, 3). 

Having constructed the couplings between the mediator and the quarks in the gauge basis, 
one still has to transform to the mass eigenstate basis. In the case of Q the final result of this 
transformation is obvious, because the Lagrangian is simply the quark part of the Yukawa sector 
of the SM. One finds 

T ^ ^ ^ (di T Vi didj^ , (2) 

where h is the physical Higgs field and = y/2m^lv with v ~ 246 GeV the vacuum expectation 
value (VEV) of H that breaks the electroweak symmetry. The lesson to learn from this exercise 
is that in order to construct MFV simplified models that describe s-channel exchange of spin-0 
resonances, the portal couplings to the SM fermions should be of Yukawa type. The above line of 
reasoning will be applied to the simplihed models in Section 


2.1.2. MFV Spin-1 s-Channel Models 

The second example that we want to discuss is even simpler than the first one. We consider 
the interactions of DM in form of a Dirac fermion x with the SM quarks through the exchange 
of spin-1 mediators which we call Z'. MFV does not restrict the couplings between the mediator 
and DM, and as a consequence the interactions take the generic form Z'^ PlX~^9^X1^PrIk) 

with Pl,r = ( 1 T 75)/2 denoting left-handed and right-handed chiral projectors. Since the bilinears 
qX^Q: and are all flavor singlets, we do not have to invoke the Yukawa couplings Y“ 

and Y'^, and simply write 


T D ^ \^gl {ui'j'^PLUi + dix'^Pidi) + g\ua^PRUi + g^da^PRdi 
i 


( 3 ) 


In fact, this expression holds both in the gauge as well as the mass eigenstate basis as long as 
the coefficients ^1^, 5 ^, and gR are flavor independent. Notice that ^ contains the case of pure 
vector or axialvector quark couplings as a special case, i.e. 9l = 9r = 9r or g^ = —gR = —Pr, 
respectively. Spin-1 s-channel simplified models of MFV type will be discussed in Section 


2.1.3. Comment on Non-MFV Models 

For the sake of argument let us also consider an example of a simplified model that does 
not conform with MFV. As a toy-model we take a Z' boson that couples vectorial to the quark 
gauge eigenstates, but differently to the first, compared to the second and third generations. We 
parameterize this non-universality by a real parameter Ay, and restrict ourselves to down-type 
quarks writing 

£ D ^ {gv -V Aydii) diX^di. 
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( 4 ) 



To go to the mass eigenstate basis requires rotating the left-handed and right-handed quark fields 
by 3 X 3 unitarity matrices These rotations will leave the term proportional to gy flavor 

diagonal, but will induce flavor off-diagonal terms of the form 


L D ^ (Lijdi'y^Pidj + Rijdi'j^PRdj) . 


( 5 ) 


with 

L = C/f diag ( 1 , 0 , 0 ) i? = C/^tdiag(l, 0 , 0 )t/|. ( 6 ) 

At this point we have to make some assumptions about the flavor structure of the ultravio¬ 
let (UV) complete model that gives rise to Q to progress further. Since the right-handed rotations 
are not observable in the SM, we assume that is the 3x3 unit matrix I 3 . This implies 
that R = diag (1,0,0) and thus there are no FCNCs in the right-handed down-quark sector. In 
contrast, the left-handed rotations are observable in the SM, because they combine to give the 
CKM matrix, i.e. V = A possible simple choice that satisfies this requirement is = V 

and U1 = I 3 , resulting in 

/iKdp v:,Vus v:,Vub\ 

L=\V*,Vud IK.P , (7) 

v:,Vus IKfelV 

which implies flavor violation in the down-type quark sector. Note that choosing Uf^ = I 3 and 
~ would give L = diag (1, 0, 0). However, this choice does not solve the new-physics “flavor 
problem”, because it is easy to see that FCNCs would then appear in the up-type quark sector. 

Using 0 it is straightforward to calculate the FCNCs induced by tree-level Z'-boson exchange. 
For instance, the new-physics amplitude relevant for kaon mixing can be estimated to be 

A{,d ^ Z' ^ U) ~ ^ 5 X 10-^ ^ . ( 8 ) 

with the mass of the Z' boson. This result should be compared to the dominant SM contribution 
to K~K mixing, which arises from top-VF boxes and is given by 


A{sd —)• box —)■ sd) 


a; 


iy,Wi 


ts 




5 X 10 


-13 


256M^ 




( 9 ) 


with = g^ j (dvr) the weak coupling constant, M\y the IF-boson mass, and yt — ^ the top Yukawa. 
A rough bound on the amount of additional flavor violation /S.y /Mz' can now be obtained by simply 
requiring that Q should be smaller in magnitude than (§. It follows that 


Ay My/ 
Mz' 


< 3 X lO-*^ 


( 10 ) 


which implies that for Ay ~ 1 the Z'-boson mass Mz' should be larger than around 3 x 10^ TeV, 
because otherwise one would be in conflict with the experimental bounds on kaon mixing. In view 
of this result it should be clear that in order to allow for interesting LHC phenomenology, one has 
to require that the simplified model is MFV. In our toy model Q this is simply achieved by setting 
Ay = 0. We finally add that for very light mediators important constraints on simplified models 
can however still arise from quark-flavor physics even if their interactions are MFV (see |31) for a 
recent comprehensive discussion). 









2.2. Note about Spins 

In many cases, there will be variations of the simplified model under consideration where the 
DM is a real or complex scalar, Dirac or Majorana fermion, or even a neutral vector. In some 
cases, even simple changes such as considering a Majorana instead of a Dirac fermion can lead 
to big changes in the phenomenology of direct detection experiments and/or annihilation. The 
classical examples are that for Majorana fermions the vector coupling vanishes identically and that 
such DM particles cannot have an electric or magnetic dipole moment. In the context of simple cut- 
and-count analyses at the LHC, the precise nature of the DM particle is generically less relevant 
in the sense that it will to first order only affect the total production cross sections. Angular 
observables that are sensitive to the structure of the dark sector have however been constructed 
and studied [SaiMlIMj, but such analyses necessarily involve topologies beyond 2 —)• $rp + l. 


3. Scalar s-channel Mediator 

A scalar particle mediator can be a very simple addition to the SM. If it is chosen as a gauge 
singlet, it can have tree-level interactions with a singlet DM particle that is either a Dirac or 
Majorana fermion, or DM that is itself a scalar. The spin-0 mediator could still be chosen as either 
a real or complex scalar, which are distinguished by the fact that a complex scalar contains both 
scalar and pseudoscalar particles, whereas the real option contains only the scalar field. We will 
consider here two choices for DM simplified models: one where the interaction with the SM is 
mediated by the real scalar, and the second where we consider only a light pseudoscalar (assuming 
that the associated scalar is decoupled from the low-energy spectrum). 

Couplings to the SM fermions can be arranged by mixing with the SM Higgs. Such models have 
intriguing connections with Higgs physics, and can be viewed as generalizations of the Higgs portal 
to DM. The impact on Higgs physics is discussed in Section |4.2| below. The most general scalar 
mediator models will of course have renormalizable interactions between the SM Higgs and the 
new scalar cj) or pseudoscalar a, as well as 4>/a interactions with electroweak gauge bosons. Such 
interactions are model-dependent, often subject to constraints from electroweak precision tests, 
and would suggest specialized searches which cannot be generalized to a broad class of models 
(unlike, for instance, the -|- j searches). As a result, for this class of simplified models with 
spin-0 mediators, we suggest to focus primarily on the couplings to fermions and the loop-induced 
couplings to gluons. The possibility that the couplings to the electroweak sector can also lead to 
interesting DM phenomenology should however be kept in mind, and can be studied in the context 
of Higgs portal DM. 


3.1. Fermionic DM 

MFV dictates that the coupling of a scalar to the SM fermions will be proportional to the 
fermion masses. However, it allows these couplings to be scaled by separate factors for the up-type 
quarks, down-type quarks, and the charged leptons. Assuming that DM is a Dirac fermion X) 
which couples to the SM only through a scalar cj) or pseudoscalar a, the most general tree-level 
Lagrangians compatible with the MFV assumption are |23( l35] : 


C-ie 


D -gx4>XX - [guyfuiUi + gdvfdidi + geyUi^i ) , 


la 


•^fermion,a 2) igx^X'l5X rpr / ^ [jluyi -\~ guyi di'J^di -|- g(.y,i^i'yb^i j ■ 

'' i 
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( 11 ) 

( 12 ) 


f 

Here the sums run over the three SM families and we are using Yukawa couplings normalized 
as y{ = y/2m{ /v with v the Higgs VEV. We parametrize the DM-mediator coupling by , rather 
than by a Yukawa coupling y^ = \/2my-jv, since the the DM particle x most likely receives its mass 
from other (unknown) mechanisms, rather than electroweak symmetry breaking. 

The most general Lagrangians including new scalars or pseudoscalars will have a potential 
containing interactions with the SM Higgs field h. As stated above, we choose to take a more 
minimal set of possible interactions, and leave the discussions of the couplings in the Higgs sector 
to the section on Higgs portal DM. Given this simplification, the minimal set of parameters under 
consideration is 


|m^, Tn^ja, gxi 9ui Odi 9i} ■ 


(13) 


The simplest choice of couplings is gu = gd = 9 £, which is realized in singlet scalar extensions of 


the SM (see Section 4.2). If one extends the SM Higgs sector to a two-Higgs-doublet model, one 


can obtain other coupling patterns such as gu oc cot /3 and g^ oc g^ oc tan /? with tan /3 denoting the 
ratio of VEVs of the two Higgs doublets. The case gu 9 d 7 ^ 9 £ requires additional scalars, whose 
masses could be rather heavy. Eor simplicity, we will use universal couplings gv = 9 u = 9 d = 9e 
in the remainder of this section, though one should bear in mind that finding ways to test this 
assumption experimentally would be very useful. 

The signal strength in DM pair production does not only depend on the masses and 
and the couplings gt, but also on the total decay width of the mediator (/>/a. In the minimal model 
as specified by 0 and (|12[), the widths for the mediators are given by: 
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(14) 

(15) 


with 


U{r) = 


1 -|- (1 — r) arctan^ 


- ly. 


/a(r) = rarctan^ 


Vt-1 


(16) 


The first term in each width corresponds to the decay into SM fermions (the sum runs over all 
kinematically accessible fermions, Nc = 3 for quarks and Nc = 1 for leptons). The second term 
is the decay into DM (assuming that this decay is kinematically allowed). The factor of two 
between the decay into SM fermions and into DM is a result of our choice of normalization of the 
Yukawa couplings. The last term corresponds to decay into gluons. Since we have assumed that 
9 v = 9 u = 9 d = 9 i: we have included in the partial decay widths r((/)/o —)■ gg) only the contributions 
stemming from top loops, which provide the by far largest corrections given that yt 3> yt etc. At 
the loop level the mediators can decay not only to gluons but also to pairs of photons and other 
final states if these are kinematically accessible. The decay rates T{(j)/a —>■ gg) are however always 
larger than the other loop-induced partial widths, and in consequence the total decay widths T^/^ 
are well approximated by the corresponding sum of the individual partial decay widths involving 
DM, fermion or gluon pairs. Notice finally that if and gu ^ gx the total widths of p/a 

will typically be dominated by the partial widths to top quarks. 
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Figure 3: Left: Two examples of one-loop diagrams with an exchange of a c^/o mediator that provide the dominant 
contribution to a mono-jet signature. Right: A tree-level graph that leads to a -f ti signal. 


3.1.1. LHC Searches 

Under the assumption of MFV, supplemented by gv = du = dd = Qi, the most relevant couplings 
between DM and the SM arising from (11) and (12) are those that involve top quarks. Two main 
strategies have been exploited to search for scalar and pseudoscalar interactions of this type using 
LHC data. The first possibility consists in looking for a mono-jet plus missing energy signal 
where the mediators that pair produce DM are radiated from top-quark loops |36) , while the second 
possibility relies on detecting the top-quark decay products that arise from the tree-level reaction 
$rp+tt [37] . In the first paper [33] that discussed the $T~^j signal the effects of DM fermions coupled 
to heavy-quark loops were characterized in terms of effective higher-dimensional operators, i.e. with 
mediators being integrated out. The effects of dynamical scalar and pseudoscalar messengers in the 
s-channel mediating interactions between the heavy quarks in the loop and DM were computed in 
characterizing the LHC signatures for DM searches in [381 IMl IMl [Ml HI] ■ 

Final states involving top-quark pairs were considered in the articles [lH HHl HU HS |39l HI] . 
Searches for a -)- tb signal [371ll2lll5| also provide an interesting avenue to probe and 
(12), while the constraints from mono-jet searches on the scalar and pseudoscalar interactions 
involving the light quark flavors are very weak due to the strong Yukawa suppression (as discussed 
in detail in |38l I46j ). and thus are unlikely to be testable at the LHC. Scenarios where the DM-SM 
interactions proceed primarily via gluons have also been considered HTj. 

Predicting mono-jet cross sections in the simplified models © and ( [l2| ) is complicated by 
the fact that the highly energetic initial-state and/or final-state particles involved in the process 
are able to resolve the structure of the top-quark loops that generate the + j signal (see the 
left-hand side of Figure]^. Integrating out the top quark and describing the interactions by an 
effective operator of the form with the field strength tensor of QCD 

and its dual, is in such a situation a poor approximation [361 13H] • Already in 

the LHC Run I environment the —)• oo limit overestimates the exact cross sections by a factor of 
5 (40) for rriy. ~ 10 GeV (m^ ~ 1 TeV) [H] . Removing the top quark as an active degree of freedom 
becomes even less justified at 13 (14) TeV, where the and pTj selection requirements have to 
be harsher than at (7) 8 TeV to differentiate the DM signal from the SM background. In order to 
infer reliable bounds on © and ( |12| ), one therefore has to calculate the mono-jet cross section 
keeping the full top-quark mass dependence. Such calculations are now publicly available at leading 
order (LO) in MCFM [38] and at LO plus parton shower (LOPS) in the POWHEG BOX [H]. Given that 
the -|- tt{bb) signals arise in the context of (jll]) and (12) at tree level (see the right-hand side 
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of Figure]^, event generation through programs like MadGraphB |48) is possible, and UFO model 
files [39] from different groups iMiiiniiiii are available for this purpose. 

Since 0 and ( |12[ ) is a simplified DM model, it is possible that the mediator can decay into 
additional states present in the full theory that we have neglected. For example, </>/a may decay 
into new charged particles which themselves eventually decay into DM, but with additional visible 
particles that would move the event out of the selection criteria of the mono-jet or similar 
searches. Another possibility is that the mediator can also decay invisibly into other particles of the 


dark sector. In either case, the expressions for F^/^ as given in (14) and (15) are lower bounds on the 


total decay-width of the mediators. To understand how the actual value ot F^/^ influences the LHC 
sensitivity, one has to recall that for (where \/~§ is some characteristic fraction of the 

center-of-mass energy of the collider in question) and m^j^/a > 2m^, DM-pair production proceeds 
dominantly via an on-shell mediator. If the narrow width approximation (NWA) is applicable, the 
mono-jet cross section factorizes into a product of on-shell production of 4>/a times its branching 
ratio into yy, i.e. a{pp — )• = a{pp —)■ 4>/a+j) Br(4>/a —> xx)- can draw three conclusions 

from this result. First, in the parameter region where > 2m^ and F^/^ <C m^ja, a change in 
r^/a will simply lead to a rescaling of the cross section, namely a{pp —)• + j) oc l/F^/^. This 

implies in particular that kinematic distributions of simple signals will to hrst approximation 
be unaltered under variations of F^/^. Second, for parameter choices where the partial decay width 
to XX DM pairs is dominant, the cross section scales as a{pp —)• + 4 ) oc g^- If the partial decay 

width to SM particles gives the largest contribution to F^/^, one has instead a{pp —)• + oc 5 ^. 

Third, the scaling a{pp —)• + j) oc g‘^g‘1 only holds for off-shell production, which occurs for 

m^la < 2m^. Notice that for the total decay width of (^ja will have a non-trivial 

impact on the constraints that the LHC can set, since the amount of off-shell production depends 
sensitively on F^/^. 

Similarly, the total decay width effect is non-trivial when the mediator can decay into other 
particles in the invisible sector beyond the cosmologically stable DM. To apply the simplified 
models framework to these scenarios, it was proposed in |39l 130] to treat the mediator width as an 
independent parameter in the simplified model characterization. 

We now turn to the constraints on these models from non-collider experiments: thermal relic 
abundance, indirect detection, and direct detection. The hrst two resnlts can be considered to¬ 
gether, as they depend on the same set of annihilation cross sections. 


3.1.2. Thermal Cross Sections 

The thermally-average annihilation of DM through the spin-0 mediators can be calculated from 
the simplihed model 0 and ( |12| ). The resulting cross sections for annihilation into SM fermions 
are given by 


{crv){xx ^ ^ f f) = Nc 


^gigiVf 


Stt 


(m2 - 4m|)2 + m2F2 


{(^v)ixx ^ a ^ ff) = Nc 


glahjrnl 


4n [(m2 — 4m|)2 -|- m^F^] 


1 - 


1 - 



(17) 


(18) 


where T denotes the DM temperature. Notably, scalar mediators do not have a temperature- 
independent contribution to their annihilation cross section, while pseudoscalars do. As T oc v'^ 
(where v is the DM velocity), there is no velocity-independent annihilation through scalars. In the 
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Universe today u ~ 1.3 x 10 ^c, so there are no non-trivial constraints on DM annihilation from 
indirect detection in the scalar mediator model (see, however, references mm)- 

The parameter space of the pseudoscalar model, on the other hand, can be constrained by indi¬ 
rect detection. Most constraints from indirect detection are written in terms of a single annihilation 
channel, and so the constraints for the full simplified model (with multiple annihilation channels 
open) require some minor modifications of the available results. In the case at hand, good esti¬ 
mates can be obtained by considering the most massive fermion into which the DM can annihilate 
(bottom and top quarks if they are accessible), as they dominate the annihilation cross section. 
Note that, away from resonance, the total width Tq entering in (18) is relatively unimportant for 
obtaining the correct indirect detection constraints. 

The thermal relic calculation requires the same input cross sections as indirect detection. Here, 
the cross sections are summed over all kinematically available final states, and can be written as 


{(Tv) = a + bT . 


(19) 


If the DM particles are Dirac fermions, one has to include a factor of 1/2 in the averaging, because 
Dirac fermions are not their own anti-particles. In the Majorana case no such factor needs to be 
taken into account. The thermal relic abundance of DM is then 




7.88 X IQ-i^x/GeV-^ 
a 3bjxf 


( 20 ) 


where Xf = m^/Tf G [20,30] with Tf the freeze-out temperature. For reasonable early Universe 
parameters, the correct relic abundance ~ 0.11 occurs in the ballpark of 

3 X 10"2® cmVs = 2.57 x 10"^ GeV^ = a + 3b/ Xf . (21) 


Keep in mind that these equations require some modihcation when the DM-mediator system is 
on resonance. Further, recall that it is unknown whether or not DM is a thermal relic, or if the 
only annihilation process in play in the early Universe proceeds through the mediator considered in 
the simplified model. Therefore, while it is appropriate to compare the sensitivity of experimental 
results to the thermal cross section, this is not the only range of parameters of theoretical interest. 


3.1.3. Direct Detection 

In contrast to the situation discussed before, elastic scattering of DM on nucleons induced by 
(p/a exchange can be very well described in terms of an EFT. Integrating out the mediators leads 
to the expressions 


0(j) 


dxdvVq 

v/2m2 


xxqq, 


o„ 


qxdvUq 

V2m2 


xibxqibq , 


( 22 ) 


at tree level, as well as contact terms consisting of four DM or quark fields. Removing the top 
quark as an active degree of freedom generates an effective interaction between DM and gluons. At 
the one-loop level, one obtains 


Og 


(^sOxS-v 

12'Kvm?, 

9 






Og = 


^s9x9v 




(23) 


by employing the Shifman-Vainshtein-Zakharov relations [52]. At the bottom- and charm-quark 


threshold, one has to integrate out the corresponding heavy quark by again applying (23). Note 
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that this matching procedure is crucial to obtain the correct DM-nucleon scattering cross section 
associated with effective spin-0 DM-quark interactions. 

The DM scattering cross section with nuclei is then obtained by calculating the nucleon matrix 


elements of the operators (22) and (23) at a hadronic scale of the order of 1 GeV. Direct detection 


provides relevant constraints only on the scalar mediator model and not the pseudoscalar case, 
since only the operators Off, and Oq lead to a spin-independent (SI) cross section, while for Oa and 
Oq the DM-nucleon scattering turns out to be spin-dependent (SD) and momentum-suppressed. 

The scalar interactions with the nuclear targets used for direct detection are (to good ap¬ 
proximation) isospin-conserving, so that the elastic DM-nucleon cross section can be written as 
{N = n,p) 


cr^-N 




9x9v 


TT 


vm 


fN, 


(24) 




where is the DM-nucleon reduced mass = m^rriN/+ m]\f) and rriN — 0.939 GeV is 

the average nucleon mass. The form factor /^r is given by 


In = 


E 

q=u,d,s 


/: 


AT + 27 


0 . 2 , 


(25) 


where the numerical value has been obtained using ~ 0.017, ~ 0.036, ~ 0.043 [53l El] 

= 1 — J2a=u ds fh — 0.904. Notice that the constraints arising from existing and future 


and 


N 


direct limits on (24) can be evaded by assuming that x is not stable on cosmological time scales. 


but lives long enough to escape the ATLAS and GMS detectors. When comparing the bounds set 
by direct detection and the LHC, this loophole should be kept in mind. 


4. Higgs Portal DM 

DM may predominantly couple to the SM particles through the SM Higgs. There are three 
broad classes of models of this kind: 

A. The DM particle is a scalar singlet under the SM gauge group, which couples through a 
quartic interaction with the Higgs. The collider phenomenology of this DM scenario has been 
extensively studied in the literature (see for instance [551 EH EZl Ea Ea Ea EH). 

B. The DM particle is a fermion singlet under the gauge symmetries of the SM, which couples 
to a scalar boson which itself mixes with the Higgs. This model class provides a specific 
realization of the s-channel scalar mediator case discussed in Section Its implications for 
the LHC have been studied for example in [62l ES EH E5|- 

C. The DM particle itself may be a mixture of an electroweak singlet and doublet [671 ESI EH] > 

in the MSSM where it has both bino and higgsino components. Generically, this is referred 
to as “singlet-doublet” DM [70] . 

The first two cases capture important features of models [ZIlEQlES] where the SM is extended 
to be classically scale invariant [7311711 Ha [76] with the aim of addressing the electroweak gauge 
hierarchy problem. 
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4-1- Scalar Singlet DM 

In the case where an additional real scalar singlet x is the DM candidate, the Lagrangian of 
the scalar Higgs portal can be written as 


•^scalar,H D ^pX. 


(26) 


where H denotes the usual SM Higgs doublet. Augmenting the Lagrangian with a discrete Z 2 
symmetry that takes x ~X H ^ H leads to stable DM, and in addition guarantees that 
there is no singlet-Higgs mixing, which leaves the couplings of the SM Higgs unaltered at tree level. 
The self-coupling of the scalar x is in general irrelevant to determining how well the portal 
coupling Ap can be probed through LHC DM searches, and thus may be ignored. 

For rrih > the most obvious manifestation of the interactions (26) is through their contri¬ 


butions to the invisible decay of the Higgs. The corresponding decay width reads 

1/2 


xx) = 


27rm/j 


1 - 


4m^ 

ml 


(27) 


with the Higgs mass and v its VEV. In fact, both ATLAS |77j and CMS m have already 
interpreted their Run I h —invisible searches in terms of the Higgs portal scenario (26). For DM 
candidates with < 10 GeV these searches are competitive with or even stronger than the SI 
results provided by direct detection experiments. 

When mh < the Higgs cannot decay on-shell to a pair of x particles, so that DM pair 

production necessarily has to proceed off-shell. The cross section for this process is then suppressed 
by an additional factor of Ap as well as the two-body phase space, leading to a rate that rapidly 
diminishes with m^. This feature makes a LHC discovery challenging even at 14TeV and high 
luminosity |61j . 


4 . 2 . Fermion Singlet DM 

A simple model including both a real scalar mediator s and a fermion DM singlet X; which 
couple through a Higgs portal is given by 

-^fermion,// 7) //gS XgS 2/xXX® ApS \H\ , (28) 

where is a Yukawa coupling in the dark sector, while the iXp and Ap terms provide the Higgs 
portal between the dark and the SM sectors. The precise values of the Higgs potential parameters 
yLs and A^ do not play an important role in the DM phenomenology at the LHC and therefore all 
features relevant for our discussion can be captured within the restricted framework yg = \s = 


In general the Higgs potential in (28) develops nontrivial VEVs for both H and s, but in order 
to keep the expressions simple it is assumed in the following that (s) = 0. The main physics 
implications are unaffected by this assumption. As a result of the portal coupling yp, the Higgs 
and the real scalar fields mix, giving rise to the physical mass eigenstates hi and /12: 


cos 9 sin 9 
— sin 9 cos 9 


(29) 


The mixing angle is defined such that in the limit 9 
Analytically, one has 

tan(20) = 


0 the dark sector is decoupled from the SM. 


2vyp 


+ XpV^ - mi 


m‘ 


(30) 
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Figure 4: Left: One-loop diagrams with an exchange of a hi/h 2 mediator that induces a mono-VF and mono-^ 
signal. Right: Two possible one-loop graphs that contribute to a mono-Higgs signal. 


while the masses of hi and /12 are given by rrih^ ~ rrih and ~ {ml + The state hi can 

therefore be identified with the bosonic resonance discovered at the LHC. 

To make contact with the scalar mediator model described in Sectionwe consider the Yukawa 
terms that follow from (28). After electroweak symmetry breaking and rotation to the mass eigen¬ 
state basis, one finds 


1 


L D-^ (cos 6 hi — sin 

v2 


0 ^ 2 )'^ y/ff - (sin 9hi+ cos 9 /i 2 ) VxXX 
f 


(31) 


Identifying /i 2 with the field (j) in & one sees that as far as the couplings between /12 and the SM 
fermions are concerned, the interactions (31) resemble those of the scalar mediator model described 


in Section with Qu = dd = 9e = 9v = — sin0. The coupling between DM and the mediator, called 
9x in 0. is instead given by = Vx cos 9. 

Another important feature of (31) is that the effective Yukawa coupling between hi and the 


SM fermions is not yj but yf cos 9. In fact, the universal suppression factor cos 9 appears not only 
in the fermion couplings but also the hiW~^W~ and hiZZ tree-level vertices as well as the loop- 


induced higg, /ii 77 , and hijZ couplings. The mixing angle and hence (28) is therefore subject to 


the constraints that arise from the ATLAS and CMS measurements of the signal strengths in Higgs 
production and decay. Global fits [79l|80] to the LHC Run I data find sin6l < 0.4. Constraints on 
the mixing angle also derive from the oblique parameters T {aka the p parameter) and S [63], but 
they are typically weaker than those that follow from Higgs physics. 


Like in the case of the scalar singlet DM model discussed before, the model (28) allows for invis¬ 
ible decays of the Higgs boson, if this is kinematically possible, i.e. mh^ > 2m^. The corresponding 
decay rate is 

2 \ 3/2 


r(/ii ^ xx) = 


2 • 2 
yt sm^ I 




Stt 


1 - 


4mz 
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(32) 


hi 


After the replacements sin 9 —>■ cos 9 and mh-^ —)• m /12 the same expression holds in the case of h 2 , 


if it is sufficiently heavy. In order to determine from (32) the invisible Higgs branching ratio, one 
has to keep in mind that all partial widths of hi to SM particles are suppressed by cos^ 9 and that 
depending on the mass spectrum also the decay hi —?■ /i 2 /i 2 may be allowed. 

Turning our attention to the signals, an important observation to make is that the phe¬ 
nomenology of the fermion singlet DM scenario is generically richer than that of the scalar mediator 
model ©• First of all, since the Lagrangian ( [28| ) leads to couplings between the scalars hi and /12 
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to massive gauge bosons as well as DM pairs, mono-ID and mono-Z signals will arise at tree level. 
The relevant diagrams are shown on the left-hand side in Figure The resulting amplitudes for 
mono-VF or mono-Z production at the LHC take the following schematic form 


A{pp —> $T + W/Z) oc sin(20) 


Sxx - ^h2 + r h2 


(33) 


where denotes the invariant mass of the DM pair and T/i^ and are the total decay widths 
of the scalars. Note that the contributions from virtual / 11 /L 2 exchange have opposite sign in (33). 


This implies that the + WjZ signal cross sections can depend sensitively on and rriy^ as 
well as the cuts imposed in the analysis. The destructive interference between the contributions of 
the two scalar mediators is also at work for mono-jets and it is well-known |621 [63l [M] that it can 
be phenomenologically relevant in direct detection. 

A second interesting consequence of (28) is that this Lagrangian gives rise to a mono-Higgs 
signal |65L I66j . Two examples of Feynman graphs that provide a contribution are given on the 
right in Figure]^ Notice that while a $x + ^ signal can also arise in the simplified s-channel scalar 
mediator scenario discussed in Section the presence of the two scalar states hi and L 2 and the 
existence of trilinear Higgs vertices such as are likely to change the mono-Higgs phenomenology 
of (28) compared to ([TT). 


4-3. Singlet-Doublet DM 

Singlet-doublet DM scenarios are the simplest example of models where the interactions between 
DM and the SM arise from mixing of a singlet with electroweak multiplets. A fermion singlet y and a 
pair of fermion doublets with opposite hypercharge denoted by = (V'l, and ^^2 = {i’t 
are introduced. Assuming that the new fields are odd under a Z 2 symmetry under which the SM 
fields are even, the Lagrangian reads 

>Cfermion,SD = i {x$X + + '02T^V’2) “ ]^msX^ “ 'niD'4^1'4^2 

- yixHiIJi - ?/ 2 X^^V 2 + b.c., (34) 


where denotes the covariant derivative. The model generalizes the bino-higgsino sector of the 
MSSM in the decoupling limit. In fact, the Yukawa couplings yi and 2/2 are free parameters, whereas 
in the MSSM they are related to the U{1)y gauge coupling. 

After electroweak symmetry breaking, singlet and doublets mix. The physical spectrum consists 
of a pair of charged particles (x"'~,x~) with mass mn and three neutral eigenstates defined by 
(xi,X2jX 3)^ = C/(x, V'lJ'02)'^! where U is the unitary matrix that diagonalizes the mass matrix 
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ms 

yiv 

V2 

y2V \ 

V2 



yiv 

Y2 

0 

mo 

(35) 

V 

y2V 

Y2 

mn 

« ) 



The DM candidate is the lightest eigenstate xij whose composition in terms of gauge eigenstates 
is xi = t^iiX + Ui2'4^i + f^i 3 V’ 2 ' singlet-doublet scenario, DM couples to the Higgs boson h 

and the SM gauge bosons through its doublet components. The induced interactions can be read 
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off from 


£ D - hxi{cl^^^.PL + ChxiXjPR)Xj - Z^Xil^{czxix,PL - ^ZxiXjPR^Xj 

- ^{U^^W-XirPLX-^ - Ut 2 W-xa^PRX+ + h.c.), (36) 


where i,j = 1,3 and 

czxiXi = - Ui2U*2 ), Chx.xj = ^{yiUi2Uji + y2UaUji ), (37) 


with g the SU{2)l coupling and cos 9w the cosine of the weak mixing angle. Due to these in¬ 
teractions, DM can annihilate to SM fermions via s-channel Higgs or Z-boson exchange and to 
bosons again through a Higgs or a Z boson in the s-channel or via Xi or irr the t-channel. 
Likewise, Higgs (Z-boson) exchange leads to SI (SD) DM nucleon scattering. The corresponding 
phenomenology has been studied in ISZIEEIEOIEI]. 

As in the case of the other Higgs portal models, a possible collider signature is the invisible 
width of the Higgs, if the decay h —>• xiXi is kinematically allowed: 


r(/i ^ XiXi) 


ruh 

dvr 




3/2 


\chx 


ixil 


(38) 


Since the Z boson couples directly to pairs of DM particles xi, the model (34) will also give rise 
to an additional contribution to the invisible decay width of the Z boson of the form 


r(-Z' xixi) 


mz A _ 4rr4 
Ovr I 


3/2 

l^ZxiXll ’ 


(39) 


if mz > ^RT-xi ■ This possibility is constrained by the Z-pole measurements performed at LEP 
which require r(Z — )• xiXi) ^ 3MeV. 

Since the model ( |34[ ) contains one charged and two neutral fermions in addition to the DM state 
Xi, LHC searches for electroweak Drell-Yan production allow to set bounds on the new fermions 
arising in scalar-doublet scenarios. The relevant production modes are qq —>■ XiXj and qq —)■ x'''X~ 
via a Z boson or qq^'^ —)• x^Xi through IT-boson exchange. Generically, the latter production mode 
has the most relevant LHC constraints. Production in gluon-gluon fusion gg XiXi through an 
intermediate Higgs produced via a top-quark loop is also possible. Like in the case of electroweakino 
production in the MSSM, final states involving leptons and provide the cleanest way to probe 
singlet-doublet models [671 [68l [81] . A particularly promising channel is for instance pp —)■ X^X 2,3 —^ 
VL^Xi-^Xi that leads to both a and 3£ + signature. The scenario (34) predicts further 

collider signals with such as mono-jets that await explorations. 


5. Vector s-Channel Mediator 

5.1. Model-Building Aspects 

One of the simplest ways to add a new mediator to the SM is by extending its gauge symmetry 
by a new (7(1)', which is spontaneously broken such that the mediator obtains a mass My [831184j . 
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Depending on whether DM is a Dirac fermion x or a complex scalar tp, the interactions this new 
spin-1 mediator take the form [T8l [85l |2T1 |86l [871 EH] 

Aermiony D V^xi^{g^ - g^-f5)x + - 9fl5)f , (40) 

>Cscaiary D ig^Vfj,{(p* ^ Vf,fj>^{gJ - gfx5)f , (41) 

f=qAi' 

where q,i and v denote all quarks, charged leptons and neutrinos, respectively. Under the MFV 
assumption the couplings of V to the SM fermions will be flavor independent, but they can depend 
on chirality (such that g^ ^ 0). For Majorana DM, the vector coupling g^ vanishes, while a real 
scalar cannot have any CP-conserving interactions with V. 

In the literature, one often finds a distinction between vector mediators with vanishing axi- 
alvector couplings {g^ = 0) and axialvector mediators with vanishing vector couplings {g^ = 0). 
Neglecting the couplings to neutrinos, the relevant parameters in the former case are 

My, g^, gl, g^, gj) , (42) 

while, in the latter case, the corresponding set is 

{m^. My, g^, g^, gj, gf] . (43) 

Note, however, that it is rather difficult to engineer purely axialvector couplings to all quarks, while 
being consistent with the SM Yukawa interactions and MFV (as explained below). In the following, 
we will consider the general case with non-zero vector and axialvector couplings. Although in this 
case the spin-1 mediator is not a parity eigenstate, we will refer to it as a vector mediator for 
simplicity. 

5.1.1. The Higgs Sector 

The most straightforward way to generate the mass of the vector mediator is by introducing an 
additional dark Higgs field $ with a non-zero VEV. Generically, this particle will not couple directly 
to SM fermions, but it could in principle mix with the SM Higgs, leading to a phenomenology similar 
to that of Higgs portal models described in Section The mass of the dark Higgs cannot be very 
much heavier than that of the vector mediator, and so may need to be included in the description 
if My is small compared to the typical energies of the collider. 

Moreover, if the theory is chiral, i.e. if g^ ^ 0, the dark Higgs will also be responsible for 
generating the DM mass. In order for the Yukawa interaction ^XX to b® gauge invariant, we have 
to require that the U (1)' charge of the left-handed and the right-handed component of the DM field 
differ by exactly Ql — Qr = Q^- Consequently, the axialvector coupling of DM to the mediator will 
necessarily be proportional to g#. The longitudinal component of V (i.e. the would-be Goldstone 
mode) then couples to x with a coupling strength proportional to g^m^/Mv. Requiring this 
interaction to remain perturbative gives the bound 

. \/47r , ^ 

implying that the DM mass cannot be raised arbitrarily compared to the mediator mass. 
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(44) 



A similar consideration also applies in the visible sector. If the axialvector couplings to the SM 
states are non-zero, the only way to have SM Yukawa couplings is if the SM Higgs doublet H 
carries a charge qn under the new gauge group. This charge must satisfy g'qu = —gu = 9d ~ 9e 
(where g' is the gauge coupling of the 17(1)0 hr order for quark and charged lepton masses to be 
consistent with the 17(1)' symmetry. However, having qn 0 generically implies corrections to 
electroweak precision measurements, so that one must require My ^ 2 TeV for consistency with 
low-energy data. 


5.1.2. Mixing with SM Gauge Bosons 

As soon as there are fermions charged under both the SM gauge group and the new 17(1)0 loop 
effects will induce mixing between the new vector mediator and the neutral SM gauge bosons, in 
particular kinetic mixing of the form 

/:kinetic 3 I , (45) 


where 7))^, = d^Vu — and = d^B^, — denote the 17(1)' and 17(l)y field strength 

tensors. Parametrically, this mixing is given by 
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<? 


{9q 


A\2 


IGvr^ 
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(46) 


If My is too close to the Z-boson mass Mz, this mixing can lead to conflicts with electroweak 
precision observables [89l EU 1901 El]- For example, the correction to the p parameter, Ap = 
M^/{M‘^cos^ 6w) — 1, can be estimated to be 


Ap ~ 


Ml 

Ml-Ml' 


(47) 


Requiring Ap < 10 ^ [92] then implies < 1 and My > 100 GeV. 


5.2. Phenomenological Aspects 

The first observation is that in models with s-channels mediators, the possibility for such parti¬ 
cles to decay back to the SM is unavoidably present. This can show up as di-jets [Sn] or di-leptons 
at the LHC. Indeed the leptonic couplings gj and gf are very tightly constrained by searches for 
di-lepton resonances IH31EH]. If the quark couplings of the mediator are equally small, it becomes 
very difficult to have sizable interactions between the SM and DM and there would typically be no 
observable DM signals. We therefore focus on the case where the quark couplings of the vector me¬ 
diator are much larger than the lepton couplings, for example because the SM quarks are charged 
under the new t/(l)' while couplings to leptons only arise at loop-level (a so-called leptophobic Z' 
boson). 

For such a set-up to be theoretically consistent we must require additional fermions charged un¬ 
der the 17(1)' and the SM gauge group to cancel anomalies. The masses of these additional fermions 
are expected to be roughly of the order of My, so they can often be neglected in phenomenology, 
unless the mass of the vector mediator is taken to be small compared to the typical energy scales 
of the collider. Indeed, it is possible to construct anomaly-free models with no direct couplings 
to leptons (for example in the context of a baryonic Z' boson [931 194|h In this case, the leptonic 
couplings will not give a relevant contribution to the DM phenomenology of the model and one can 
simply set g\ = gf = 1). 
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5.2.1. Collider Searches 

If the vector mediator is kinematically accessible at the LHC, the resulting phenomenology 
depends crucially on its decay pattern. For arbitrary vector and axialvector couplings, one finds in 
the case of Dirac DM the following expression for the total width: 


by 


My 

I2n 


i=f,X 


4m? \ 


{gY? + {gtf + 


A\2 


2 

m,- 

K 




(48) 


Here the sum extends over all fermions i that are above threshold, while iV* = 3 for quarks and 
= 1 for leptons and DM. 

There are several important conclusions that can be drawn from (48). The first one concerns the 
maximal size that the couplings can take to be consistent with Ty/My < 1, which is a necessary 
requirement in order for a perturbative description of the mediator to be valid. Assuming that 


My ^ rrii and setting for simplicity Qq = = g and g^ = gf = 0, one finds that Ty/My ~ D.5g^. 

This implies that one has to have o' < 1.4 in order for the width of the mediator to be smaller than 


its mass and values of g significantly below unity for the NWA (which calls for Ty/My < 0.25) to 
be applicable. 

In cases where the NWA can be used, production and decay factorize such that for instance 
a{pp — )■ Z + xx) = (y{pp -A Z + V) X Br(I4 —)• xx)- The resulting LHC phenomenology is thus 
determined to first approximation by the leading decay mode of the vector mediator. Considering 
a situation with My ::?> m* and gj = g^ = {), one finds that decays into quarks dominate if 
3x1 gq ~ while invisible decays dominate if ^/g^ > 4. For g^/g^ — 4 both decay channels 
have comparable branching ratios. If invisible decays dominate, the strongest collider constraints 
are expected from searches for in association with SM particles. To illustrate this case, we 
discuss mono-jet searches below. If, on the other hand, the invisible branching ratio is small, we 
expect most of the mediators produced at the LHC to decay back into SM particles. In this case, 
strong constraints can be expected from searches for heavy resonances, and we focus on di-jet 
resonances. 


Mono-Jets 

LHC searches for $rp plus jet signals place strong constraints on the interactions between quarks 
and DM mediated by a vector mediator [2011211 na [Ml sa [951 Ezi 198] . The corresponding cross 
sections can be calculated at next-to-leading order (NLO) in MCFM |38j and at NLO plus parton 
shower in the POWHEG BOX |46j . Some of the relevant diagrams are shown in Figure If the 
mediator is too heavy to be produced on-shell at the LHC and assuming equal vector couplings 
of the mediator to all quarks as well as gj = gf = 0, the mono-jet cross section at the LHC 
is proportional to ( 5 ^)^ (S'x same scaling applies if the mediator is forced to be off-shell 

because My < 2m^ so that decays into DM are kinematically forbidden. 

For <C My <C \/s, with ^/s the center-of-mass of the collider, the mediator can be pro¬ 
duced on-shell and subsequently decay into a pair of DM particles. If the mediator width is small 
enough for the NWA to be valid, the mono-jet cross section will be proportional to the product 
{g ^)^ Br(F —)• xx). If we fix the ratio g^/g^, the invisible branching ratio will be independent of an 
overall rescaling of the couplings, so that we simply obtain (T{pp —)• -T j) oc {gYY- If we rescale 
only one of the couplings, on the other hand, the resulting change in the mono-jet cross section 
will depend on the dominant decay channels of the mediator. If the total width of the mediator is 
dominated by its decays into quarks, the mono-jet cross section will be invariant under a rescaling 
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Figure 5: Left: An example of a LO diagram that leads to mono-jet events through the s-channel exchange of a 
spin-1 vector resonance V. Right: At the NLO level both virtual and real corrections have to taken into account in 
order to obtain a infrared finite result. 


of the quark coupling , because the change in the production cross section is compensated by the 
change in the invisible branching ratio. If, on the other hand, invisible decays dominate, both the 
production cross section and the invisible branching ratio will be invariant under a (small) change 
in the coupling g^. 

The same general considerations apply for axialvector couplings instead of vector couplings. 
In particular, the production cross section of the vector mediator is largely invariant under the 
exchange g^ <->■ g^. Note, however, that for —>• My/2 the phase space suppression is stronger 

for axialvector couplings than for vector couplings, such that for ~ My/2 the monojet cross 
section is somewhat suppressed for a mediator with purely axialvector couplings. 

In many situations invisible decays and decays into quarks will both lead to a non-negligible 
contribution to Ty as given in (48) and furthermore this width may become so large that one 
cannot use the NWA to derive simple scaling laws. If becomes close to My/2 there can also 
be contributions from both on-shell and off-shell mediators. As a result, all relevant parameters 
{ruy., My, g^ and g^) must in general be taken into account in order to calculate mono-jet cross 
sections. 


Di-Jets 

Searches for di-jet resonances exploit the fact that any mediator produced from quarks in the 
initial state can also decay back into quarks, which lead to observable features in the distribution 
of the di-jet invariant mass and their angular correlations. However, for small mediator masses 
the QCD background resulting from processes involving gluons in the initial state completely over¬ 
whelms the signal. The most recent di-jet searches at the LHC therefore focus mostly on the region 
with di-jet invariant mass rrijj > 1 TeV. For smaller mediator masses, the strongest bounds are in 
fact obtained from searches for di-jet resonances at UA2 and the Tevatron [98]. An interesting op¬ 
portunity to make progress with the LHC even in the low-mass region is to consider the production 
of di-jet resonances in association with other SM particles, such as W or Z bosons, which suffer 
from a significantly smaller QCD background |99L llOOj . 

An important complication concerning searches for di-jet resonances results from the fact that 
the width of the mediator can be fairly large. The steeply falling parton distribution functions then 
imply that the resonance will likely be produced at lower masses, leading to a significant distortion 
of the expected distribution of invariant masses rrijj. Existing searches for narrow resonances 
therefore typically do not apply to vector mediators with couplings of order unity. Nevertheless, 
the shape of the resonance can still be distinguished from SM backgrounds and it is still possible 
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to constrain these models using specifically designed searches |98] . 

A number of such searches have been considered in |98j . The central conclusion is that, at 
least for < 1, bounds on My become stronger as is increased, because the enhancement 
of the production cross section is larger than the reduction of the detection efficiency resulting 
from the increasing width. Indeed, there are still stringent bounds on mediators with width as 
large as Ty ~ My/ 2. It is crucial to take these bounds into account when interpreting DM 
searches at the LHC in terms of simplihed models with an s-channel vector mediator, because 
they apply to a wide range of models and in many cases complement or even surpass other search 
strategies. A promising strategy to constrain even broader resonances may be to study di-jet angular 
correlations, such as the ones considered in the context of constraining four-fermion operators (see 
for instance ununnn). 


5.2.2. Direct Detection 

Depending on the coupling structure of the vector mediator, the interactions of DM with nu¬ 
clei can proceed via SI or SD scattering off nucleons. The corresponding cross sections at zero- 
momentum transfer are given by 


^si 

^X-N 


/y2 

^^X-N f2 


^SD 

^X-N 


irMy 
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ajv 


(49) 


where N stands for either p or n, while /jv and ajy denote the effective nucleon couplings, 
take the form 


fp — 5x i‘^9u + 9d) ^ 


fn — fl'x (9u + ‘^9d ) ) 


They 

(50) 


and 


ap,n = 9x gf . (51) 

q=u,d,s 

The coefficients encode the contributions of the light quarks to the nucleon spin. They are 

given by [92] 

Au^P") = aS'^^ = 0.84 ± 0.02 , 

= -0.43 ± 0.02 , (52) 

As^P'^ = As^^'^ = -0.09 ± 0.02 . 


Potential cross terms such as are suppressed in the non-relativistic limit (either by the mo- 

mentum transfer or the DM velocity, both of which lead to a suppression of 10“^ or more), and 
can therefore be neglected. 

Substituting the expressions for the effective couplings into the formulas for the DM-nucleon 
scattering cross sections, we obtain 
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Crucially, SI interactions receive a coherent enhancement proportional to the square of the target 
nucleus mass, leading to very strong constraints from direct detection experiments unless the DM 
mass is very small. Consequently, the estimates above imply that for Qq — 1, SI interactions are 
sensitive to mediator masses of up to My ~ 30 TeV, while SD interactions only probe mediator 
masses up to around My ~ 700 GeV. This should be contrasted with the constraints arising from 
the LHC, which are close to identical for vector and axialvector mediators. 


5.2.3. Annihilation 

Two processes contribute to DM annihilation in the early Universe: annihilation of DM into SM 
fermions and (provided My < direct annihilation into pairs of mediators, which subsequently 
decay into SM states. For the first process, the annihilation cross section is given by 
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where Ty is the total decay width of the vector mediator as given in (48). For ~ My/2 the 


annihilation rate receives a resonant enhancement, leading to a very efficient depletion of DM. 

An important observation is that for = 0, the annihilation cross section is helicity-suppressed. 
For nih <C < mt the factor rn/^/m^ can be very small, such that it is important to also include 
the p-wave contribution for calculating the DM relic abundance. Including terms up to second 
order in the DM velocity v, we obtain for the special case = 0 the expression 


icrv){xX ^ qq) = 
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Finally, the annihilation cross section for direct annihilation into pairs of mediators is given by 


{(^v){xx VV) 



6{gi)H9^f + {9^)^]M^y 


(57) 


We note that for the coupling strengths and mass ranges typically considered in the context of 
LHC DM searches, it is easily possible to achieve sufficiently large annihilation cross sections to 
deplete the DM abundance in the early Universe. In fact, the generic prediction in large regions of 
parameter space would be that the DM particle is underproduced. In this case, the observed DM 
relic abundance can still be reproduced if one assumes an initial particle-antiparticle asymmetry 
in the dark sector, such that only the symmetric component annihilates away and the final DM 
abundance is set by the initial asymmetry. 
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6. t-Channel Flavored Mediator 


If the DM is a fermion Xi the mediator can be a colored scalar or a vector particle 4>. We focus 
on the scalar case, which makes contact with the MSSM and is easier to embed into a UV-complete 
theory. A coupling of the form (pxQ requires either y or (/> to carry a flavor index in order to 
be consistent with MFV. We choose the case where the colored scalar cj) carries the flavor index 
(much like in the MSSM case, where the colored scalar quarks come in the same flavors as the SM 
quarks). This class of models has been considered previously in |104l 1106111031110511107111081 [T6] , 
while models where x carries the flavor index have been studied in [TIMfTTnifTTT] . 

There are variations where the mediator couples to right-handed up-type quarks, right-handed 
down-type quarks, or left-handed quark doublets. For definiteness, we discuss the right-handed up- 
type case (the other cases are obtained in a similar fashion). In this case, there are three mediators 
(j)i = {ft, c, t}, which couple to the SM and DM via the interaction 

■^fermion,ti —^ E + h.c. (58) 

i=l,2,3 


Note that MFV requires both the masses Mi^ 2,3 of the three mediators to be equal and universal 
couplings g = 51 , 2,3 between the mediators and their corresponding quarks Ui = {u,c,t}. This 
universality can however be broken by allowing for corrections to (58) and the mediator masses 
which involve a single insertion of the flavor spurion Because of the large top-quark Yukawa 

coupling, in this way the mass of the third mediator and its coupling can be split from the other 
two. In practice this means that the generic parameter space is five-dimensional: 


{m^, Ml, 2 , Ms, 51,2, 53 } • (59) 

These simplified models are very similar to the existing ones for squark searches [112j . and results 
can often be translated from one to the other with relatively little work. Note that most studies 
will involve 51,2 together with Mi,2 or 53 together with M 3 . So specific applications will often have 
a smaller dimensional space of relevant parameters. In the discussion below, we restrict attention 
to the parameter space with 51,2, Mi,2, and rriy’. For models where 53 and M3 are relevant, 
see [rmfTTilfTTTlfl^ . 


6.1. Collider Constraints 

Given the masses and couplings, the widths of the mediators are calculable. One finds 
r{4,i ^ xui) = 

X Mf -b -b — 2M?m^, — 2Mfm‘^ — 

Mi,m^ > niui . 

Mi > m^,mui . 


(1 - 

4;M,, 


Mf 


(60) 


Unless the final-state quark Ui is a top quark, the given limiting cases are always very good ap¬ 
proximations to the exact widths. 

In the context of (58) the production channels that lead to a + 3 signal are uu —)> XX + 9^ 
ug XX + s-iid ug —)> XX + u- Examples of the relevant Feynman diagrams are shown on the 
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Figure 6: A + j signal can arise in the t-channel mediator scenario from initial-state gluon emission (left) and 
associated mediator production (middle). Initial-state gluon splitting processes and gluon emission from the t-channel 
mediator is also possible but not shown. Pair production of the mediator u in gluon fusion leads instead to $rp -\- 2j 
events (right). Quark-fusion pair production either via s-channel gluon or t-channel DM exchange also contributes 
to the latter signal. 


left and in the middle of Figure In addition, if the colored mediator u is sufficiently light it 
may be pair produced from both gg or uu initial states. This gives rise to a signature 

as illustrated by the graph on the right-hand side in the same figure. If the DM particle is a 
Majorana fermion also the uu and uu initial states contribute to the production of mediator pairs. 
The latter corrections vanish if x is a Dirac fermion. From this brief discussion, it should be clear 
that t-channel models can be effectively probed through both mono-jet and squark searches. 

6.1.1. Mono-Jet Searches 

Given that in all recent mono-jet analyses a second hard jet is allowed, the corresponding 
LHC searches are sensitive in t-channel models to the contributions not only from initial-state 
gluon radiation and associated production, but also to mediator pair production. Since the relative 
importance of the different channels depends on , Mi, and gi as well as the imposed experimental 
cuts, all corrections should be included in an actual analysis. General statements about the leading 
partonic channel are however possible. For what concerns events the diagram in the middle 

of Figure usually gives the dominant contribution. Compared to uu —)• XX + 9-, this process 
benefits from a phase-space enhancement, the larger gluon luminosity, and the fact that jets from 
initial state radiation tend to be softer. If the mass Mi is small, diagrams with gluon emission from 
the mediator can also be important, but these graphs are subdominant if the mediator is heavy, 
since they are 1/M^ suppressed. Notice that the dominance of the associated production channel 
is a distinct feature of f-channel models that is not present in the case of s-channel mediators, nor 
is it relevant in supersymmetric theories where the mediator is a squark. The relative importances 
of the different -\- j and + 2j channels depend sensitively on how gi compares to the strong 
coupling constant gs. In the limit gi <C gs, pure QCD pair production dominates, while in the 
opposite case graphs with DM exchange are more important. Detailed studies of the bounds on the 
coupling gi as a function of Mi and that arise from Run I mono-jet data have been presented 

in [imfTn^ . 

6.1.2. Squark Searches 

If the t-channel mediator is light it can be copiously produced in pairs at the LHG and then 
decay into DM and a quark. The resulting phenomenology is very similar to squark pair production 
in the MSSM with a decoupled gluino. There is however one important difference which has 
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to do with the fact that in supersymmetric theories the coupling between the squarks and the 
neutralino x is necessarily weak. The cross section for squark pair production through t-channel 
exchange of DM is therefore negligible. This is not the case in t-channel mediator scenarios, because 
gi is a free parameter and thus it is possible to enhance significantly the u pair production rate 
associated to t-channel DM exchange. As already mentioned, the quark-fusion pair production 
cross section depends on whether x is a Dirac or a Majorana particle. In the former case only 
uti-initiated production is non-zero, while in the latter case also the uu and uu initial states furnish 
a contribution. The constraints from LHC squark searches on t-channel mediator models with both 
Dirac and Majorana DM have been investigated thoroughly in |106l 1105111071 [TOH] . These studies 
show that squark and mono-jet searches provide comparable and complementary bounds on a wide 
range of the parameter space of f-channel scenarios, depending on the masses of the mediator and 
DM. Especially in the case where the DM particle and the mediator are quasi-degenerate in mass, 
mono-jet searches turn out to be superior. 


6.2. Scattering with Nucleons 

Away from resonance and neglecting light quark-mass effects, the SI scattering cross section of 


Dirac DM and nucleons that is induced by (58) reads 
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Here fp = 2 and fn = ^ and hence the SI cross sections for protons and neutrons are different in 
the t-channel scenario. Using the same approximations the subleading SD scattering cross section 
takes the form 
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with the numerical values for given in (52). Notice that for Majorana DM, the SI scattering 

cross section vanishes and the expression for is simply obtained from ( |62[ ) by multiplying the 

above result by a factor of 4. 

Since in t-channel models with Dirac DM one has 0, the existing direct detection 

constraints dominate over the collider bounds up to very low DM masses of around 5 GeV. For 
Majorana DM instead — as a result of the lack of the enhancement from coherence in DM-nucleus 
scattering — the LHC constraints turn out to be stronger than the direct detection limits for DM 
masses up to of a few hundred GeV. 


6.3. Annihilation Rates 

The main annihilation channel of DM in the framework of (58) is XX UiUi. In the Dirac case 
this leads to a s-wave contribution of the form 
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if quark masses are neglected (remember that an additional factor of 1 /2 has to be included in the 
thermal averaging). In the Majorana case, annihilation to SM quarks is instead p-wave suppressed 
and given for zero quark masses by 


(o-u)(xx um) = 
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In the parameter space where the mediator and the DM particle x a-re quasi-degenerate 
in mass and the ratio (Mj — is comparable to or below the freeze-out temperature, co¬ 

annihilation effects become important mm- For both Dirac and Majorana fermions the 
annihilation cross section for x4’i '*^*5 can be written as 
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if quark-mass effects and u^-suppressed terms are neglected. In addition the mediators (j)i can 
self-annihilate. While for both Dirac and Majorana DM annihilation to gluons 
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proceeds via s-wave, the process —>■ UiUi is p-wave suppressed and hence subdominant. Finally, 

for Majorana DM the reaction —>■ UiUi (and its charge conjugate) is possible. The relevant 

s-wave contribution in this case reads 


{av){<j)i<j)i UiUi) = 
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Assuming that the relic abundance is thermally produced, one finds that for Dirac DM 

there is no region in the parameter space that satisfies the combined constraints arising from the 
LHC searches, direct detection, and Q.yh?'. Therefore the simple model (58) with Dirac DM cannot 


be regarded as a complete model in describing the interactions between the dark and the visible 
sectors. In the case of Majorana fermions satisfying all three requirements is possible, but the mass 
of DM must be larger than about 100 GeV. If DM is lighter there must be other channels for DM 
to annihilate into, which calls for additional new physics. 


7. Conclusions 

The primary goal of this document is to outline a set of simplified models of DM and their 
interactions with the SM. It can thus serve as a summary and proposal for the simplified models 
to be implemented in future searches for DM at the LHC. The list of models discussed includes 
spin-0 and spin-1 s-channel mediator scenarios as well as t-channel models. The most important 
prototypes of Higgs-portal scenarios are also described. To motivate our choice of simplihed models, 
a number of guiding principles have been given that theories of DM-SM interactions should satisfy 
in order to be useful at LHC energies (and possibly beyond). Based on these criteria building further 
simplified (or even complete) DM models is possible. While the focus is on giving a brief account 
of the LHC signals that seem most relevant in each of the simplified models, we have also provided 
expressions and formulas for reference that allow the reader to derive the constraints from direct 
and indirect searches for DM. There is still useful work to be done to improve our understanding 
of simplified DM models, and room to devise creative new searches that can discover or constrain 
them. 

While most of the discussion in this work centers around simplified models, we emphasize that 
the given examples represent only “theoretical sketches” of DM-SM interactions, and that they 
by no means exhaust the whole spectrum of possibilities that the DM theory space has to offer. 
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They are neither meant to form self-contained, complete pictures of DM interactions at the LHC, 
nor are they meant to be model-independent and general enough to cover the entirety of the DM 
landscape. In order to do justice to the range of options in the DM theory space, it is thus important 
when searching for DM at the LHC to frame the results of searches in terms of all three types of 
theoretical frameworks: EFTs, simplified models, and UV complete theories. Only in this way is 
it possible to maximize the search coverage for DM at LHC Run H, and have the largest possible 
impact on our understanding of the particle properties of DM. Simplified models thus play a crucial 
role in this endeavor. 
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